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ABSTRACT 

We studied all blazars of known redshift detected by the Fermi satellite during its first three- 
months survey. For the majority of them, pointed Swift observations ensures a good multi- 
wavelength coverage, enabling us to to reliably construct their spectral energy distributions 
(SED). We model the SEDs using a one-zone leptonic model and study the distributions of 
the derived interesting physical parameters as a function of the observed 7-ray luminosity. 
We confirm previous findings concerning the relation of the physical parameters with source 
luminosity which are at the origin of the blazar sequence. The SEDs allow to estimate the 
luminosity of the accretion disk for the majority of broad emitting line blazars, while for the 
line-less BL Lac objects in the sample upper limits can be derived. We find a positive corre- 
lation between the jet power and the luminosity of the accretion disk in broad line blazars. In 
these objects we argue that the jet must be proton-dominated, and that the total jet power is 
of the same order of (or slightly larger than) the disk luminosity. We discuss two alternative 
scenarios to explain this result. 

Key words: BL Lacertae objects: general — quasars: general — radiation mechanisms: non- 
thermal — gamma-rays: theory — X-rays: general 



1 INTRODUCTION 

The Large Area Telescope (LAT) onboard the Fermi Satellite, in 
the first months of operation succeeded to double the number of 
know 7-ray emitting blazars (Abdo et al. 2009a, hereafter A09). 
Excluding radio-galaxies, in the list of the LAT 3-months survey 
we have 104 blazars, divided into 58 Flat Spectrum Radio Quasars 
(FSRQs, including one Narrow Line Seyfert 1 galaxy), 42 BL Lacs 
and 4 sources of uncertain classification. Of these, 89 have a known 
redshift, and a large fraction of sources have been observed with the 
Swift satellite within the 3 months survey, while others have been 
observed by Swift at other epochs. 

The combination of Fermi and Swift data are of crucial impor- 
tance for the modelling of these blazars, even if the LAT data are an 
average over the 3 months of the survey, so that, strictly speaking, 
we cannot have a really simultaneous optical to X-ray SED, even 
when the Swift observations have been performed during the three 
months of the survey. Despite this, the optical/UV, X-ray and 7-ray 
coverage can often define or strongly constrain the position and the 
flux levels of both peaks of the non-thermal emission of blazars. 
Furthermore, as discussed in Ghisellini, Tavecchio & Ghirlanda 
(2009, hereafter Paper 1), the data from the UV Optical Telescope 
(UVOT) of Swift are often crucial to disentangle the non-thermal 
beamed jet emission from thermal radiation produced by the accre- 
tion disk. When the accretion disk is visible it is then possible to 
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estimate both the mass of the black hole and the accretion luminos- 
ity and to compare it with the power carried by the jet. In Paper I 
we did this study for the 23 most luminous blazars, exceeding, in 
the 7-ray range, an observed luminosity of 10*** erg s^^. In Tavec- 
chio et al. (2009, hereafter T09) we considered the BL Lac objects 
detected by Fermi with the aim to find, among those, the best can- 
didates to be TeV emitters. Here we extend these previous studies 
to the entire Fermi blazar sample of the 3 months survey. Besides 
finding the intrinsic physical properties characterising the emitting 
region of these Fermi blazars, like magnetic field, particle density, 
size, beaming factor and so on, the main goal of our study is to in- 
vestigate if there is a relation between the jet power and the disk 
luminosity, as a function of the 7-ray luminosity. 

This implies two steps. First we collect the data to 
construct the SED of all sources, taking advantage of pos- 
sible Swift observations, that we analyse, and of archival 
data taken from the NASA Extragalactic Database (NED, 
'http : / /nedwww . ipac . caltech . edu/ Secondly, we ap- 
ply a model to fit these data, that returns the value of the physi- 
cal parameters of the emitting region, and, if possible, the value of 
the accretion disk luminosity and of the black hole mass. We do 
this for all the 85 sources (out of 89) with a minimum coverage 
of the SED. Among these blazars, there are several BL Lac ob- 
jects for which there is no sign of thermal emission produced by 
a disk, nor of emission lines produced by the photo-ionising disk 
photons. For them we derive an upper limit of the luminosity of a 
standard, Shakura-Sunyaev (1973) disk. We also derive the distri- 
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butions of the needed physical parameters, and their dependence 
on the 7-ray luminosity and on the presence/absence of a promi- 
nent accretion disk. As expected, the jets of line-less BL Lacs have 
less power than broad line FSRQs, have bluer spectra, and their 
emitting electrons suffer less radiative cooling. This confirms ear- 
lier results (Ghisellini et al. 1998) explaining the so called blazar 
sequence (Fossati et al. 1998). 

The main result of our analysis concerns the relation between 
the jet power and the accretion disk luminosity. We find that they 
correlate in FSRQs, for which we can estimate both the black hole 
mass and the accretion luminosity. We discuss two alternative sce- 
narios to explain this behaviour. For BL Lacs (with only an upper 
limit on the accretion luminosity), we suggest that the absence of 
any sign of thermal emission, coupled to the presence of a relatively 
important jet, strongly suggests a radiatively inefficient accretion 
regime. 

In this paper we use a cosmology with h = 57a = 0.7 and 
JIm ~ 0.3, and use the notation Q — Qx 10^ in cgs units (except 
for the black hole masses, measured in solar mass units). 



2 THE SAMPLE 

Tab. [T] lists the 89 blazars with redshift of the A09 catalogue of 
Fermi detected blazars. Besides their name, we report the redshift, 
the K-corrected 7-ray luminosity in the FennifLKT band, (see e.g. 
Ghisellini, Maraschi & Tavecchio 2009), if there are observations 
by the Swift satellite, if the source was detected by EGRET and the 
classification. For the latter item we have maintained the same clas- 
sification as in A09, but we have indicated those sources classified 
as BL Lac objects but that do have broad emission lines. 

The distinction between BL Lacs and FSRQs is traditionally 
based on the line equivalent width (EW) being smaller or larger 
than SA. With this definition, we classify as BL Lacs those sources 
having genuinely very weak or absent lines and also objects with 
strong lines but whose non-thermal continuum is so enhanced to 
reduce the line EW. This second category of "BL Lacs" should 
physically be associated to FSRQs. We have tentatively made this 
distinction on Tab. [T] based on some information from the litera- 
ture about the presence, in these objects, of strong broad lines. In 
the rest of the paper we will consider them as FSRQs. 

The source PMN 0948-0022 belongs to still another class, be- 
ing aNarrow Line Seyfert 1 (NLSyl; FWHM~1500 km s"^). This 
source has been discussed in detail in Abdo et al. (2009b) and Fos- 
chini et al. (2009): its SED and general properties are indistinguish- 
able from FSRQs and we then assign this source to this group. 

We then have 57 FSRQs, 1 NLSyl and 31 BL Lac objects. Of 
the latter, 6 have strong broad lines or have been also classified as 
FSRQs (by Jauncey et al. 1989). Four of the FSRQs (indicated with 
the superscript n in Table [T) do not have a sufficient data coverage 
to allow a meaningful interpretation, and we will not discuss them. 
The 23 blazars (21 FSRQs and 2 BL Lacs, but with broad emission 
lines) with an average 7-ray luminosity exceeding lO**** erg s^^ 
written in italics have been discussed in Paper I. The 25 underlined 
BL Lacs ai-e discussed in T09. One of these, (08I4-l-425=OJ 425) is 
shown also here, since it can be fitted both with a pure synchrotron 
self-Compton (SSC) and an "External Compton" (EC) model (see 
next section). 

In summary, in this paper we present the SED and the cor- 
responding models for 37 blazars: 32 FSRQs and 5 BL Lacs that 
are suspected to be FSRQs with broad lines hidden by the beamed 
continuum. However, when discussing the general properties of the 
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Table 1. The 89 Fermi blazars in the A09 catalogue with redshift, including 
two blazai-s that have no z in the A09 lists (0332^03 and 1553+1 1). In the 
last 4 columns we indicate the logarithm of the average 7-ray luminosity as 
observed by Fermi during the first 3 months of survey (cgs units); if there 
are Swift observations; if the source was detected by EGRET (UL stands for 
an upper limit given by Fichtel et al. 1994); the classification of the source 
(B=BL Lac; Q=FSRQs; NL=Narrow line Seyfert galaxy; U=uncertain clas- 
sification), a: Narrow Line Seyfert 1, analysed in Abdo et al. (2009b) and 
Foschini et al. (2009). b: redshift uncertain, n: not studied in this paper 
due to lack of multiwavelength data; *; defined as BL Lacs for the EW of 
the lines, but broad emission lines are present. **; classified as FSRQ in 
Jauncey et al. (1989). 

Fermi blazars, we will consider the entire blazar sample of Tab.[T] 
with the only exception of the 4 FSRQs with a very poor data cov- 
erage. These are 85 sources. To this aim we will use the results of 
Paper 1, concerning the most 7-ray luminous blazars, and we will 
apply our model also to the BL Lacs shown in T09. For the ease of 
the reader, we report in Tableland Table[5]the physical parameters 
of all the 85 blazars. 
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Table 1. - continue - 



3 SWIFT OBSERVATIONS AND ANALYSIS 

For 33 of the 37 blazars studied in this paper there are Swift obser- 
vations, with several of them being observed during the 3 months 
of the Fermi survey. The data were analysed with the most re- 
cent software SWIFT_REL3 . 2 released as part of the Heasof t 
V. 6.6.2. The calibration database is that updated to April 10, 
2009. The XRT data were processed with the standard procedures 
(XRTPIPELINE V. 0.12. 2). We considered photon counting 
(PC) mode data with the standard 0-12 grade selection. Source 
events were extracted in a circular region of aperture ~ 47", 
and background was estimated in a same sized circular region far 
from the source. Ancillary response files were created through the 



xrtmkarf task. The channels with energies below 0.2 keV and 
above 10 keV were excluded from the fit and the spectra were re- 
binned in energy so to have at least 30 counts per bin. Each spec- 
trum was analysed through XSPEC with an absorbed power-law 
with a fixed Galactic column density from Kalberla et al. (2005). 
The computed errors represent the 90% confidence interval on the 
spectral parameters. Tab. [2] reports the log of the observations and 
the results of the fitting the X-ray data with a simple power law 
model. 

UVOT (Roming et al. 2005) source counts were extracted 
from a circular region 5"— sized centred on the source position, 
while the background was extracted from a larger circular nearby 
source-free region. Data were integrated with the uvotimsum 
task and then analysed by using the uvotsource task. The ob- 
served magnitudes have been dereddened according to the formu- 
lae by Cardelli et al. (1989) and converted into fluxes by using stan- 
dard formulae and zero points from Poole et al. (2008). Tab.[3]list 
the observed magnitudes in the 6 filters of UVOT. 



4 THE MODEL 

We use the model described in detail in Ghisellini & Tavecchio 
(2009, hereafter GT09). It is a relatively simple, one zone, homoge- 
neous synchrotron and Inverse Compton model, aiming at account- 
ing the several contributions to the radiation energy density pro- 
duced externally to the jet, and their dependence upon the distance 
of the emitting blob to the black hole. Besides the synchrotron ra- 
diation produced internally to the jet, we in fact consider radiation 
coming directly from the disk (i.e. Dermer & Schlickeiser 1993), 
the broad line region (BLR; e.g. Sikora, Begelman & Rees 1994), 
a dusty torus (see Blazejowski et al. 2000; Sikora et al. 2002), the 
host galaxy light and the cosmic background radiation. 

The emitting region is assumed spherical, of size rdiss, moving 
with a bulk Lorentz factor F and is located at a distance i?diss from 
the black hole of mass Al. The bolometric luminosity of the accre- 
tion disk is Ld. The jet accelerates in its inner parts with F oc R^^^ 
(R is the distance from the black hole), up to a value Fmax- In the 
acceleration region the jet is parabolic (following, e.g. Vlahakis & 
Konigl 2004) and beyond this point the jet becomes conical with a 
semi-aperture angle if) (assumed to be 0. 1 for all sources). 

The energy particle distribution A'^(7) [cm~"^] is calculated 
solving the continuity equation where particle injection, radiative 
cooling and pair production (via the 7-7 process) are taken 

into account. The injection function (5(7) [cm^"^ s^^] is assumed 
to be a smoothly joining broken power-law, with a slope Q{-~f) oc 
7~°i and 7^*2 below and above a break energy 7b: 



Q(7) 



(7/7b)" 



1 + (7/7b)-''i+' 



(1) 



The total power injected into the source in the form of rel- 
ativistic electrons is P( — nicC^V J Q{'y)'ydj, where V = 
(47r/3)rdij,j, is the volume of the emitting region. 

The injection process lasts for a light crossing time rdiss/c, 
and we calculate N{-y) at this time. This assumption comes from 
the fact that even if injection lasted longer, adiabatic losses caused 
by the expansion of the source (which is travelling while emitting) 
and the corresponding decrease of the magnetic field would make 
the observed flux to decrease. Therefore our calculated spectra cor- 
respond to the maximum of a flaring episode. 

The BLR is assumed for simplicity to be a thin spherical 
shell located at a distance -Rblr = 10^^ L\^^^ cm. A fraction 
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^■'-1.0 


0.2/3 


2.0+0.4 


0.29+0.05 


0332^03C') 


25/02/2009 


1.38 


1.4±0.2 


2/6 


3.5+0.4 


2.4+0.2 


0537^41 


08/10/2008 


3.94 


1.7±0.1 


23/20 


8.24+0.9 


4.63+0.47 


0650+453* 


14/02/2009 


8.64 


2.1 ±0.6 


5/5 


0.76+0.11 


0.28+0.04 


0716+332 


27/01/2008 


5.93 


1.8±0.3 


3/3 


0.9+0.2 


0.4+0.2 


0948+0022** 


05/12/2008 


5.22 


1.8it0.2 


5/8 


4.4+0.2 


2.2+0.1 


0954+556 


05/03/2009 


0.853 


0+11 

"•^-1.8 


0.1/2 


2.2+0.5 


1.9+0.5 


1030+61* 


03/07/2009 


0.63 


1.9±0.7 


0.3/4 


0.62+0.10 


0.27+0.04 


1055+018* 


19/07/2008 


4.02 


1.8±0.5 


1/3 


3.5 


1.8 


1057-79 


20/01/2008 


8.76 


1.9+0.15 


6/10 


2.8+0.4 


1.3+0.3 


1127-145 


24/03/2007 


4.04 


1.31+0.05 


116/94 


9.8+0.5 


7.2+0.4 


1144-379* 


21/11/2008 


7.5 


1.96+0.32 


1/5 


1.0+0.3 


0.46+1.2 


1156+295 


21/11/2008 


1.68 


1.52+0.14 


11/7 


2.7+0.4 


1.7+0.2 


1226+023 


10/05/2008 


1.79 


1.56+0.05 


88/68 


833+19 


521+26 


1244-255 


17/01/2007 


6.85 


1.8+0.5 


3/4 


2.1+0.5 


1.0+0.4 


1253-055('') 


08/08/2008 


2.12 


1.76+0.04 


88/94 


10.5+0.2 


5.5+0.1 


1308+32 


20/08/2008 


1.27 


1.6+0.2 


5/10 


3.0+0.2 


1.8+0.1 


1508-055f ) 


20/02/2007 


6.09 


1.8+0.2 


8/6 


0.95+0.08 


0.50+0.04 


1510-089 


10/01/2009 


7.8 


1.4+0.11 


194/316 


6.4+1.3 


4.4+0.4 


1803+784 


17/02/2007 


4.12 


1.5+0.1 


14/14 


2.8+0.4 


1.8+0.2 


1846+322C') 


28/12/2008 


9.93 


1.7+0.2 


10/7 


0.78+0.05 


0.45+0.03 


1849+67('=) 


16/11/2006 


4.66 


1.5+0.1 


22/14 


2.6+0.1 


1.59+0.07 


1 908-20 K-') 


08/03/2007 


9.24 


1.4+0.1 


5/10 


2.7+0.1 


1.83+0.05 


1920-211 


07/08/2007 


5.69 


1.6+0.3 


1/3 


1.4+0.1 


0.87+0.09 


2141+175(9) 


19/04/2007 


7.35 


1.71+0.06 


94/67 


2.23+0.04 


1.47+0.03 


2144+092('') 


28/04/2009 


4.56 


1.6+0.2 


13/9 


1.6+0.1 


1.01+0.07 


2155+31C) 


01/04/2009 


7.42 


1.0+0.6 


0.6/1 


1.0+0.1 


0.8+0.1 


2201+171(') 


08/12/2009 


4.56 


i.8otni 


13/17 


1.22+0.05 


0.62+0.03 


2204-54 


21/12/2008 


1.72 


1.6+0.1 


9/17 


2.2+0.1 


1.37+0.07 


2230+1 14(™) 


19/05/2005 


4.76 


1.50+0.05 


74/66 


6.5+0.1 


4.2+0.1 


2345-155r) 


23/12/2008 


1.64 


1.6+0.4 


0.3/1 


0.39+0.05 


0.23+0.03 



Table 2. Results of the X-ray analysis. *: poorly determined spectrum, the C-Statistic was used to fit the spectrum. Therefore, the x2/ dof indicates instead the 
C-stat value and the PHA bins. ** Abdo et al. (2()09b). a: Average of 2 observations the same day. h: Average of 8 observations between 8 and 20 August 2008. 
c: Average of 4 observations between 26/02/2007 and 13/09/2008. d: Average of 2 observations on 28/12/2008 and 20/02/2009. e: Average of 3 observations 
between 16/11/2006 and 04/08/2008. /■ Average of 4 observations between 8 and 17 March 2007. g: Average of 9 observations between 19/04/2007 and 
15/01/2009. h: Average of 2 observations the same day. i: Average of 3 observations between 1 and 17 April 2009. /; Average of 7 observations between 
08/12/2006 and 18/04/2009. m: Average of 5 observations between 19/05/2005 and 01/05/2007. n: Average of 2 observations on 23/12/2008 and 10/01/2009. 



/blr = 0.1 of the disk luminosity is re-emitted by broad lines. 

1/2 

Since JJblr oc , the radiation energy density of the broad Une 
emission within the BLR is constant, but is seen amplified by a fac- 
tor ~ by the moving blob, as long as iidiss < ^blr- A dusty 
torus, located at a distance Rm = 2.5 x 10^^ L^^^ cm, reprocesses 
a fraction /ir (of the order of 0.1-0.3) of La through dust emission 
in the far IR. Above and below the accretion disk, in its inner parts, 
there is an X-ray emitting corona of luminosity Lx (we almost al- 
ways fixes it at a level of 30% of L^). Its spectrum is a power law 
of energy index ax = 1 ending with a exponential cut at Ec =150 
keV. The specific energy density (i.e. as a function of frequency) of 
all these external components are calculated in the comoving frame, 
and used to properly calculate the resulting External inverse Comp- 
ton (EC) spectrum. The internally produced synchrotron emission 
is used to calculate the synchrotron self Compton (SSC) flux. 



5 SOME GUIDELINES FOR THE MODELLING 

In this section we follow (and somewhat repeat) the arguments pre- 
sented in Paper 1, adding some considerations for line-less BL Lac 
objects. 

A general comment concerns the flux at low (sub-mm to ra- 
dio) frequencies. The one-zone homogeneous model here adopted 
is aimed to explain the bulk of the emission, and necessarily re- 
quires a compact source, self-absorbed (for synchrotron) at ~ 10^^ 
Hz. The flux at radio frequencies must be produced further out 
in the jet. Radio data, therefore, are not directly constraining the 
model. Indirectly though, they can suggest a sort of continuity be- 
tween the level of the radio emission and what the model predicts 
at higher frequencies. 



5.1 Strong line objects 

Consider first sources whose inverse Compton flux is dominated by 
the EC process with photons of the broad line region. 
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source 


OBS date 


V 


B 


U 


Wl 


M2 


W2 


0133+47 


18/11/2008 


15.89±0.02 


1646±0.01 


15.89±0.01 


16.30±0.02 


16.56±0.03 


16.74±0.02 


0208-512 


29/12/2008 


17.69±0.07 


18.06±0.05 


17.01±0.04 


16.79±0.03 


16.69±0.03 


17.03±0.03 


0218+35 


12/12/2008 


>194 


>20.3 


>19.9 


>20.3 


>20.1 


>21.0 


0332-403 


25/02/2009 


17.00±0.07 


17.59±0.05 


16.85±0.05 


17.23±0.06 


17.29±0.07 


18.20±0.07 


0537^1 


08/10/2008 


16.04±0.02 


1648±0.02 


15.78±0.02 


16.01±0.02 


16.02±0.02 


16.26±0.02 


0650+453 


14/02/2009 


>19.8 


20.1 ±0.2 


19.3±0.2 


19.3±0.1 


19.4±0.1 


19.8±0.1 


0716+332 


27/01/2008 


— 


- 


- 


— 


- 


17.71±0.03 


0948+0022* 


05/12/2008 


18.2±0.2 


18.56±0.07 


17.79±0.06 


17.48±0.06 


17.50d=0.06 


17.55d=0.05 


0954+556 


05/03/2009 


17.6±0.2 


17.97±0.06 


16.91±0.04 


16.80±0.04 


16.8it0.2 


16.94±0.08 


1030+61 


03/07/2009 


184±0.1 


19.14±0.08 


18.36±0.07 


18.56±0.07 


18.5±0.1 


19.08±0.08 


1055+018 


19/07/2008 


— 


— 


— 


— 


16.97±0.04 


— 


1057-79 


20/01/2008 


— 


— 


— 


— 


17.31±0.03 


17.36±0.01 


1127-145 


24/03/2007 


1646±0.02 


16.70±0.01 


15.64±0.01 


15.51±0.01 


15.55±0.01 


15.79±0.01 


1144-379 


21/11/2008 


— 


— 


— 


— 


— 


20.0±0.1 


1 156+295 


21/11/2008 


— 


- 


- 


16.85±0.01 


- 


- 


1226+023 


10/05/2008 


12.72±0.01 


12.92±0.01 


11.92±0.01 


11.48±0.01 


11.33±0.01 


11.34±0.01 


1244-255 


17/01/2007 


— 


- 


— 


- 


17.44±0.03 


— 


1253-055 


08/08/2008 


16.80±0.09 


16.25±0.06 


16.41±0.05 


16.44±0.05 


16.35±0.06 


16.50±0.04 


1308+32 


20/08/2008 


— 


— 


— 


— 


16.71±0.03 


— 


1508-055 


20/02/2007 


174±0.2 


17.5±0.1 


16.56±0.06 


16.38±0.09 


16.8±0.2 


16.97±0.09 


1510-089 


10/01/2009 


16.76±0.08 


17.04±0.04 


16.18±0.02 


16.51±0.02 


1642±0.04 


16.52±0.03 


1803+784 


17/02/2007 


16.51±0.03 


16.99±0.02 


16.32±0.02 


16.50±0.02 


16.55±0.02 


16.65±0.01 


1846+322 


28/12/2008 


18.9±0.1 


19.27±0.09 


18.50±0.07 


1841±0.06 


1845±0.06 


18.62±0.05 


1849+67 


16/11/2006 


17.6±0.1 


17.93±0.06 


17.28±0.06 


17.36±0.06 


16.13±0.03 


- 


1908-201 


08/03/2007 


16.8±0.1 


17.40±0.09 


16.80±0.08 


17.09±0.09 


17.34±0.08 


17.66±0.09 


1920-211 


07/08/2007 








16.73±0.03 






2141+175 


19/04/2007 


16.16±0.02 


16.39±0.02 


15.21±0.02 




15.09±0.03 


15.12±0.03 


2144+092 


28/04/2009 


18.2±0.2 


1848±0.09 


17.67±0.07 


17.50±0.06 


17.68±0.07 


18.00±0.05 


2155+31 


01/04/2009 








214it04 




>20.9 


2201+171 


08/12/2009 


16.93±0.06 


17.71±0.05 


17.23±0.04 


17.62±0.05 


17.7±0.1 


17.86±0.05 


2204-54 


21/12/2008 


17.89±0.08 


18.18±0.04 


16.98±0.03 


16.83±0.03 


16.99±0.04 


1740±0.04 


2230+114 


19/05/2005 






16.52±0.02 




16.38±0.04 


16.70±0.03 


2345-155 


23/12/2008 


18.5±0.1 


18.61±0.06 


17.91±0.05 


17.57±0.07 


17.66±0.05 


17.78±0.04 



Table 3. UVOT Observed magnitudes. * see also Abdo et al. (2009b) and Foschini et al. (2009). 



• When the UVOT data define an optical-UV bump, we inter- 
pret it as the direct emission from the accretion disc. This assump- 
tion allows us to determine both the black hole mass and the accre- 
tion rate. The maximum temperature (and hence the vFu peak of 
the disk luminosity) occurs at ~5 Schwarzschild radii and scales as 
Tmax oc {LilLEddf'^M-'^''^. The total optical-UV flux gives 
id [that of course scales as (Ld/I/Edd) A/]- Therefore we can 
derive both the black hole mass and the accretion rate. For good 
UVOT data, the method is sensitive to variations of less than a fac- 
tor 2 in both the black hole mass and the accretion rate (see the 
discussion and Fig. 2 in GhiselUni et al. 2009a). 

• The ratio between the high to low energy emission humps 
(Lc/Ls) is directly related to the ratio between the radiation 
to magnetic energy density U!^/U^. In this case the assumption 
i?BLR = 10^^Lj^45 cm gives 

where we have assumed that Ul- ^ [^blr- 

• The peak of the high energy emission (vc) is produced by the 
scattering of the line photons (mainly hydrogen Lyman-a) with 
electrons at the break of the particle distribution (7peak). Its ob- 
served frequency is fc ~ 2!/ij,ar(57peak/ (1 + A steep (energy 



spectral index a > 1) spectrum indicates a peak at energies below 
100 MeV, and this constrains TSj^g^i^. 

• Several sources whose Compton flux is dominated by the 
EC component may have rather small values of 7peak. Electrons 
with these energies may emit, by synchrotron, in the self-absorbed 
regime. In these cases the peak of the synchrotron component is the 
self-absorption frequency. 

• In powerful blazars the radiative cooling rate is almost com- 
plete, namely even low energy electrons cool in a dynamical 
timescale rdiss/c. We call 7cooi the random Lorentz factor of those 
electrons halving their energies in a timescale fdiss/c. When the 
EC process dominates, 7cooi is small (a few). Therefore the corre- 
sponding emitting particle distribution is weakly dependent of the 
low energy spectral slope, si, of the injected electron distribution. 

• The strength of the SSC relative to the EC emission depends 
on the ratio between the synchrotron over the external radiation 

energy densities, as measured in the comoving frame, (7s'/f^cxt- 
Within the BLR, Uc^t depends only on F'^, while Us depends on 
the injected power, the size of the emission, and the magnetic field. 
The larger the magnetic field, the larger the SSC component. The 
shape of the EC and SSC emission is different: besides the fact that 
the seed photon distributions are different, we have that the flux at 
a given X-ray frequency is made by electron of very different ener- 
gies, thus belonging to a different part of the electron distribution. 
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In this respect, the low frequency X-ray data of very hard X-ray 
spectra are the most constraining, since in these cases the (softer) 
SSC component must not exceed what observed. This hmits the 
magnetic field, the injected power (as measured in the comoving 
frame) and the size. Conversely, a relatively soft spectrum (but still 
rising, in vFv) indicates a SSC origin, and this constrains the com- 
bination of B, rdiss and P/ even more. 



5.2 Line-less BL Lacs 

When the SED is dominated by the SSC process, as discussed by 
Tavecchio et al. (1998), the number of observables (peak fluxes and 
frequencies of both the synchrotron and SSC components, spectral 
slopes before and after the peaks, variability timescale) is sufficient 
to fix all the model parameters. The sources whose high energy 
emission is completely dominated by the SSC process are line- 
less BL Lac objects, with less powerful and less luminous jets. The 
lack of external photons and the weaker radiation losses make 7cooi 
much larger in these sources than in FSRQs. The low energy slope 
of the injected electrons, si, in these cases coincides with the low 
energy slope of the emitting distribution of electrons. 



5.3 Upper limits to the accretion luminosity and mass 
estimate 

For line-less BL Lac objects we find an upper limit to the accretion 
disk luminosity by requiring that the emission directly produced by 
the disk and by the associated emission lines are completely hidden 
by the non-thermal continuum. By assuming a typical equivalent 
width of the emission lines observed in FSRQs (~ 100 A), and the 
one defining BL Lacs (EW< SA), we require that the disk lumi- 
nosity is a factor at least ~ 20 below the observed non-thermal 
luminosity. Note that we assume that the disk is a standard Shakura 
& Sunyaev (1973) geometrically thin optically thick disk. That this 
is probably not the case will be discussed later. 

A crude estimate of the black hole mass is obtained in the 
following way. Assume that the dissipation region, in units of 
Schwarzschild radii, is the same in BL Lacs and FSRQs. Then 
an upper limit to the mass is derived assuming that light crossing 
times do not exceed the typical variability timescales observed. A 
lower limit can be derived by considering the value of the magnetic 
field in the emitting region. Small black hole masses imply smaller 
dimensions, thus smaller Poynting flux (for a given B-field). To 
avoid implausibly small values of it, we then derive a lower limit 
to the black hole mass. Together, this two conditions give masses 
in the ranee 10* - 10^ Mq. As we discuss later (§6.1 and Tab.lU, 
for BL Lacs in which independent mass estimates are available in 
literature, the masses derived in this way agree quite well. 



6 RESULTS 

In Fig. [T] and Fig. |2] we show the distributions of the parameters 
derived from model fitting the SED of the entire sample of Fermi 
blazars (85 sources). Shaded areas in these histograms correspond 
to the 25 BL Lacs for which only an upper limit to the accretion 
luminosity could be derived (the plotted areas correspond to 26 ob- 
jects because they include two states of PKS 2155-304). They can 
be considered as "genuine" BL Lacs, namely blazars whose emis- 
sion lines are intrinsically weak or absent, and not hidden by the 
beamed continuum. 

The SED of the 37 blazar studied in this paper, together with 
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10 
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Figure 1. Distribution of the magnetic field, location of the dissipation re- 
gion in units of the Schwai'zschild radius, black hole mass, redshift, power 
Pj' injected in relativistic electrons (as measured in the comoving frame), 
bulk Lorentz factor, break energy (7b) and maximum energy (7max) of 
the injected distribution of electrons for all the 85 blazars. Shaded areas 
correspond to BL Lacs with only an upper limit on their disk accretion lu- 
minosity. 



the best fitting model, are shown in Figs. lllH20| of the Appendix. 
In Tab.|2]and Tab.[3]we list the result of our XRT and UVOT anal- 
ysis for the 33 blazars (out of the 37 blazars studied in this paper) 
with Swift data. All optical-UV fluxes shown in Fig.ll ll-|20|have 
been de-reddened according to the Galactic Av given in the NED 
database. 

Tab. |4] reports the parameters used to compute the theoretical 
SEDs and Tab. |5] lists the power carried by the jet in the form of 
radiation, electrons, magnetic field and protons (assuming one pro- 
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Table 4. Input pai'ameters used to model the SED. Sources in italics have been discussed in Paper 1. Note that -Rblr is a deiived quantity, not an independent 
input parameter It is listed for an easy comparison with -Rdiss- Col. [1]: name; Col. [2]: redshift; Col. [3]: dissipation radius in units of 10^^ cm and (in 
parenthesis) in units of /Jg; Col. [4]: black hole mass in solar masses; Col. [5]: size of the BLR in units of 10^^ cm; Col. [6]: power injected in the blob 
calculated in the comoving frame, in units of 10*^ erg s~^; Col. [7]: accretion disk luminosity in units of 10^^ erg s~^ and (in parenthesis) in units of -^Edd; 
Col. [8]: magnetic field in Gauss; Col. [9]: bulk Lorentz factor at -Rdiss; Col. [10]; viewing angle in degrees; Col. [11], [12] and [13]; minimum, break and 
maximum random Lorentz factors of the injected electrons; Col. [14] and [15]; slopes of the injected electron distribution [Q{'y)] below and above 7^,; For all 
cases the X-ray corona luminosity Lx = O.SLd- Its spectral shape is assumed to be oc exp{—hu/150 kcV). 

ton per emitting electron). For the ease of the reader, in these tables 
we report also the values found and presented in Paper 1 . 

6.1 Redshifts and black hole masses 

The redshift distribution of Fermi FSRQs extends to larger values 
than for BL Lacs (Fig.[TJ. This has been already pointed out in A09 
and is the consequence of BL Lacs being less 7-ray luminous than 
FSRQs. 

The derived black hole masses are in the range (1~ 
60)x lO^M© (Fig.[B. The object for which we could estimate the 



least massive black hole is PIVIN 0948+0022, that is a IvILSyl dis- 
cussed in Abdo et al. (2009b), where a mass of l.SxlO^M© was 
found. 

We have searched in the literature other estimates of the black 
hole masses of our blazars, and report them in Tab. [6] They have 
been mainly derived from the FWHM of the emission lines, through 
the assumption of virial velocity of the broad line clouds. There is a 
rough agreement between our and the other estimates, but note that 
for specific objects the reported estimates vary by a factor 3-10. It 
is comforting that the black hole masses assumed here for line-less 
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Table 4. - continue - 

compare this power with the power the jet carries in the form of 
bulls motion of particles and fields, since P/ is measured in the co- 
moving frame. To have comparable quantities, we should multiply 
P/byr^. 

Among FSRQs, 4 sources have JJdiss > J?blr (0215+015 
and 1520+319, discussed in Paper I, plus 0954+556 and 1622- 
253). This means a reduced radiation energy density that in turn 
implies a weaker radiative cooling and a larger 7peak. AH other 



BL Lacs are consistent with the existing estimates present in the 
literature. 



6.2 Injected power, location of the dissipation region and 
bulk Lorentz factors 

The injected power in rclativistic electrons, as measured in the co- 
moving frame, is in the range P/ = 10'*'^-10'*'* erg s~^ for FSRQs 
and a lO^^-lO"*^ erg s~^ for BL Lacs. Note that we should not 
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Figure 2. Distributions of tlie slopes of tlie injected electron distribution 
(si and S2 are the slopes before and after 71,, respectively), the value of 
the random Lorentz factor 7pcak of ths electrons radiating at the peaks of 
the SED, the value of the random Lorentz factor 7cool of electrons cooling 
in one light crossing time r^iss/c, accretion disk luminosities (also in units 
of the Eddington one), the ratio between the outflowing and accretion mass 
rate and the power spent by the jet to produce the radiation we see. Shaded 
areas coiTespond to BL Lacs with only an upper limit on their disk accretion 
luminosity. 



FSRQs dissipate within the BLR, at a distance of the order of 300- 
1000 Rs (Fig.[T}- For BL Lacs we have the same range of i?diss- 

The distribution of the bullc Lorentz factor is rather narrow, 
being contained within the 10-15 range, with few BL Lacs having 
r between 15 and 20. The average F for BL Lacs is somewhat 
larger than for FSRQs. 
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Table 5. Jet power in the form of radiation, Poynting flux, bulk motion of 
electrons and protons (assuming one proton per emitting electron). Sources 
in italics have been analysed in Paper I and their jet powers are here reported 
for completeness. 



6.3 Magnetic field 

On average, we need a slightly larger magnetic field B in FS- 
RQs (1-10 G) than in BL Lacs (0.1-1 G). Note that in our list we 
lack TeV BL Lacs not detected by Fermi (for them, see T09), that 
have more extreme properties (and smaller magnetic fields) than 
the sources considered here. Note also that there is one FSRQ with 
a small magnetic field: this is 1520+319, whose -Rdiss is at very 
large distances. 



6.4 Particle distribution 

The average properties of the injected particle distribution can be 
seen in Fig.[T]and Fig.|2] Note that the injected paiticle distribution 
(5(7) [cm~^ s"'^] is different from the emitting particle distribu- 
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Table 5. - continue - 

tion N{'y) [cm~^], that is the solution of the continuity equation 
evaluated at the light crossing time rdiss/c. It can be seen that the 
diversity of the blazar spectra requires a rather broad range of S2, 
the injected slope for the high energy electrons (between 2 and 4.5), 
steeper for BL Lac objects. Consider also that when radiative cool- 
ing is important (almost always in FSRQs, see the distribution of 
7cooi in Fig. |2}, the high energy emitting particle distribution will 
be characterised by a power law slope S2 + 1, i.e. steeper still. The 
fact that BL Lacs are characterised by a flatter 7-ray spectrum in 
the Fermi band (A09, Ghisellini, Maraschi & Tavecchio 2009) is 
not due to a flatter slope of the A'^(7) distribution above 7b, but to 
their SED peaking at energies close to or higher than a few GeV. 
Also the distribution of si is broad, with FSRQs having 



harder slopes. But in this case the radiative cooling is stronger, 
and A'"(7) cx 7"^ from 7cooi to 7b, making si less important. For 
BL Lacs, instead, the cooling is much weaker, and in several cases 
A'^(7) oc 7^^^ at intermediate energies, mostly contributing to the 
X-ray band. This requires a softer si. There are however extreme 
cases (TeV emitting BL Lacs, with an hard spectrum) present in 
T09, but not here (they are not detected by Fermi) that require a 
very hard si or even a cutoff in the electron distribution (i.e. no 
electrons injected below some critical value). 

The distributions of 7b, 7max (Fig. [Hi are clearly different for 
FSRQs and BL Lacs, with BL Lacs requiring much larger val- 
ues. This, together with the weaker cooling for BL Lacs (imply- 
ing a larger 7cooi) determines the distribution of 7pcak, the value of 
Lorentz factors of those electrons producing most of the radiation 
we see. This confirms earlier results concerning the interpretation 
of the blazar sequence (i.e. Ghisellini et al. 1998; Celotti & Ghis- 
ellini 2008). 

The relation of 7pcak with the energy densities and with 7cooi 
is shown in Fig. [3] The top panel shows 7peak as a function of 
the sum of the magnetic and radiation energy density as measured 
in the comoving frame. The grey symbols are the values for the 
blazars studied in Celotti & Ghisellini (2008). 



6.5 Disk luminosities 

The accretion disk luminosities of FSRQs (Fig |2} derived by the 
model fits are in the range 10'"^-10''^ erg s"\ and the upper limits 
for BL Lacs indicate, always, values below 10''^ erg s~^. In Ed- 
dington units, FSRQs always have values above 10~^, and BL Lacs 
always values below. Bearing in mind that for some FSRQs our es- 
timates are poor (when the beamed non-thermal flux hide the ther- 
mal emission or when Swift data are missing) this result is intrigu- 
ing. It is in perfect agreement with the "blazar's divide" between 
broad line and line-less objects proposed by Ghisellini, Maraschi 
& Tavecchio (2009). It also agrees with the scenario proposed by 
Cavaliere & D'Elia (2002) and Bottcher & Dermer (2002). We will 
further discuss this point later. 



6.6 Accretion and outflow mass rates 

Fig. [2] shows the ratio of the accretion (A/in) and outflow (Mout) 
mass rates. For FSRQs we set A/out = Pj/(rc^): therefore Afout 
suffers from the same uncertainties of Pj, derived assuming one 
proton per emitting electron and that all electrons emit. For com- 
pleteness we show this distribution also for BL Lacs, but the values 
are in this case completely uncertain, for two reasons. First, since 
only an upper limit on Ld is derived, we have a corresponding up- 
per limit for Min. On the other hand, the disks in BL Lacs may well 
be radiatively inefficient: if so, they will have larger accretion rates 
than the ones corresponding to a standard disk. We will estimate 
Afin for BL Lacs later, using the assumption of Pj — AfinC^. 

For FSRQs, instead, the distribution is more meaningful (bear- 
ing in mind the limitations mentioned above) and indicates that the 
mass outflowing rate, on average, is 1-10% of the mass accretion 
rate. The distribution is rather narrow, and this may indicate that 
the mass outflow rate of the jet (derived assuming one proton per 
emitting electron) is linked with the mass accretion rate. In other 
words, the matter of the jet may come directly from the accreting 
one, with other process, like entrainment, less important. 
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9.25 


C09 










1144-379 


8.7 


7.6 


P05 










1156+295 


9 


9.11 


C09 


8.56 


X05 


8.54 


X05 






8.63 


P05 


8.54 


L06 






1226+023 


8.9 


9.10 


F03b 


9.21 


W04 


8.6 


P05 






8.92 


L06 










1253-055 


8.9 


8.48 


W04 


7.9 


P05 


8.28 


L06 


1308+32 


8.84 


8.94 


C09 


9.24 


W04 






1502+106 


9.48 


9.50 


C09 


8.74 


L06 






1510-089 


8.84 


8.62 


W04 


8.31 


X05 










8.22 


X05 


8.20 


L06 






1551+130 


9 


9.19 


C09 










1633+382 


9.7 


9.74 


C09 


8.67 


L06 






1803+784 


8.7 


8.57 


W04 


7.92 


L06 






2141+175 


8.6 


8.98 


F03b 


8.14 


X05 


8.05 


X05 






7.95 


L06 










2227-088 


9.18 


9.40 


C09 










2230+114 


9 


8.5 


P05 










2251+158 


9 


9.10 


W04 


8.5 


P05 


8.83 


L06 



Table 6. Estimates of the mass of the black hole for the blazars in oiir 
sample. Values are given for the logarithm of the black hole mass measured 
in solar masses. We Ust here only those blazars for which we have found 
another, independent, mass estimate. The first value (Col. 2) is the estimate 
found in this paper. In the top part of the table we list BL Lacs for which 
we found only an upper limit on the disk luminosity. In this case our es- 
timates of the black hole mass is very uncertain. References: B03: Barth 
et al. (2003); C09: cheng et al. (2009); F03a: Falomo et al. (2003a); F03b: 
Falomo et al. (2003a); L06: Liu et al. (2006); P05: Plan et al. (2005); W02: 
Wu et al. (2002); W04: Wang et al. (2004); W08: Wagner (2008); Wo05: 
Woo et al. (2005). X05; Wie et al. (2005); 



6.7 Emission meclianisms 

For all but four FSRQs (0215+015, 0954+556, 1520+319 and 
1622-253) the dissipation region is within the BLR, that provides 
most of the seed photons scattered at high frequencies. The main 
emission processes are then synchrotron and thermal emission from 
the accretion disk for the low frequency parts, and External Comp- 
ton for the hard X-rays and the 7-ray part of the spectrum. The 
X-ray corona and the SSC flux marginally contribute to the soft 
X-ray part of the SED. When i?diss < Rshn the overall non- 
thermal emission is rather insensitive to the presence/absence of 
the IR torus, since the bulk of the seed photons are provided by the 
broad lines. Instead, for the 4 FSRQs with iidiss > -Rblr, the ex- 
ternal Compton process with the IR radiation of the torus is crucial 
to explain their SED. 

For BL Lacs (see the SEDs in T09) the main mechanism is 
SSC, but in few cases this process is unable to account for a very 
large separation between the synchrotron and the inverse Compton 
peaks of the SED, without invoking extremely large bulk Lorentz 
factors (larger than 100). These are the cases where the SED is 
better modelled invoking an extra source of seed photons for the 
IC process, besides the ones produced by synchrotron in the same 
zone. One possibility is offered by the spine/layer scenario (Ghis- 
ellini, Tavecchio & Chiaberge 2005), in which a slower layer sur- 
rounds the fast spine of the jet. The radiative interplay between the 
two structures enhances the IC flux and can account for the ob- 
served SED in these cases. 

Another problem, with the standard one-zone SSC scenario, 
concerns the ultrafast (i.e. minutes) variability sometimes seen at 
high energies (Aharonian et al. 2007, Albert et al. 2007). This can- 
not be accounted for by the simple models, and require other emit- 
ting zones or extra population of electrons (see e.g. Ghisellini & 
Tavecchio 2008; Ghisellini et al. 2009b; Giannios, Uzdensky & 
Begelman 2009). 

For all our sources the importance of the 7-7 process 
(which is included in our model) is very modest, and does not in- 
fluence the observed spectrum, nor the derived jet power, discussed 
below. 



7 JET POWER 

Table |5] lists the power carried by the jet in the form of radiation 
(Pi), magnetic field (Pb), electrons (Pc) and cold protons (Pp, as- 
suming one proton per emitting electron). All the powers are cal- 
culated as 

P = TrrLsr'/Scf/; (3) 

where [/,' is the energy density of the i component, as measured in 
the comoving frame. We comment below on each contribution: 

• The power carried in the form of the produced radiation, 

Pr = TrrdissF^/^c (7,!^^, can be re-written as [using U^^^ — 
L7(4^rLsC)]: 

4 4(5* 452 

where L is the total observed non-thermal luminosity (L' is in the 
comoving frame) and Ul.^^ is the radiation energy density produced 
by the jet (i.e. excluding the external components). The last equality 
assumes 6^ ^ 1 /F. This is a almost model-independent quantity, 
since it depends only on the adopted 5, that can be estimated also 
by other means, namely superluminal motions. 
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• When calculating Pe (the jet power in bulk motion of emitting 
electrons) we include their average energy, i.e. = nc{'y)mc(? ■ 
Usually, when estimating this quantity, we have the problem of de- 
termining 7min'7ieC^, the minimum energy of the electron distri- 
bution [where, for steep distribution functions N{'y), most of the 
electrons are]. Tliis problem is much alleviated here, since the as- 
sumed form of the particle injection function Q('y) is rather flat 
at low energies. The amount of the electrons at low energies then 
mainly depends on cooling, making N{'y) oc 7^^ down to 7cooi, 
and flatter below. Thus the total amount of electrons contributing to 
Pe depends on cooling, not on a pre-assigned shape of the particle 
distribution (including a pre-assigned 7min)- In this sense the Pc 
derived here is less arbitrary. Furthermore, in the case of luminous 
FSRQs, the X-ray flux can be reliably associated to the EC mech- 
anism (this occurs when the slope is very hard, because the SSC 
component tends to be rather softer, see the discussion of this point 
in Celotti & Ghisellini 2008). In this case the low energy X-ray 
data are crucial to fix 7min: a too high value makes the modelled 
X-ray flux to underestimate the observed one. Since in this sources 
the radiative cooling is severe, this agrees with the fact that 7min 
must be small, of the order of 7cooi or less. As a final comment, 
consider that the estimate of Pc includes only those electrons con- 
tributing to the emission. Since it is unlikely that all the electrons 
present in the emitting region are accelerated, Pe is a lower limit. 
On the other hand, when (7) is greater than a few, the contribution 
of the accelerated electrons to Pc may dominate over the contribu- 
tion of the cool (not accelerated) ones. 

• For Pp (the jet power in bulk motion of cold protons) we have 
assumed that there is one proton per emitting electron, i.e. electron- 



positron pairs are negligible. This is a crucial assumption. Partly, it 
is justified within the context of our model because we take into ac- 
count the pair production process, and we find that pairs are always 
negligible. If a substantial amount of pairs comes from the inner 
regions of the jet we must explain why they have survived anni- 
hilation, important in the inner, more compact and denser regions 
(Ghisellini et al. 1992). If they have survived because they were 
hot (thus they had a smaller annihilation cross section) then they 
should have produced a large amount of radiation (that we do not 
observe). In doing so, they should have cooled rapidly, and then an- 
nihilate. These considerations (see also similar comments in Celotti 
& Ghisellini 2008) lead us to accept the assumption of one proton 
per electron as the most reliable. 

For BL Lacs the presence or not of electron-positron pairs is less 
of a problem, because the mean energy of the emitting electrons is 
large, approaching the rest mass-energy of a proton. In this cases 

Pc ~ Pp. 

• Pb is derived using the magnetic field found from the model 
fitting. There can be the (somewhat contrived) possibility that the 
size of the emitting region is smaller than the one considered here, 
i.e. inside the jet there could be smaller volumes where the mag- 
netic field lines reconnect, and in this case the total Poynting flux 
of the jet can be larger than what we estimate here. We consider 
that this is unlikely. 

Fig.|4]shows the distributions of the jet powers and the bottom 
panel, for comparison, shows the distribution of the accretion disk 
luminosities. Grey shaded ares correspond to BL Lacs for which 
we could estimate only an upper limit to their disk luminosity. Fig. 
[5] shows the fraction of the total jet power (i.e. Pp 4- Pc + Pb) 
transformed in radiation (er), carried by electrons (ec) and Poynting 
flux (eb)- 

For FSRQs the power carried in radiation (Pi ) is larger than 
Pc. This is a consequence of fast cooling: electrons convert their 
energy into radiation in a time shorter than rdiss/c and the radia- 
tion component can in this time accumulate more energy than what 
remains in the electrons (even if they are continuously injected dur- 
ing this time). The distribution of Pb is at slightly smaller values 
than the distribution of Pr, indicating that the Poynting flux cannot 
be at the origin of the radiation we see. As described in Celotti & 
Ghisellini (2008), this is a direct consequence of the large values 
of the so called Compton dominance (i.e. the ratio of the Comp- 
ton to the synchrotron luminosity), since this limits the value of the 
magnetic field. 

To justify the power that the jet carries in radiation (we insist: 
it is the least controversial quantity) we are forced to consider the 
power carried by the jet in the form of protons. Following the con- 
sideration made above, the simplest and most reasonable assump- 
tion is to assume that there is one proton per electrons. If so, Pp for 
FSRQs is a factor ~ 10-100 larger than Pr, meaning an efficiency 
of 1-10% for the jet to convert its bulk kinetic motion into radia- 
tion (see also the top panel of Fig. |5]l- This is reasonable: most of 
the jet power in FSRQs goes to form and energise the large radio 
structures, and not into radiation. On the other hand, we do not have 
yet a firm handle on how much power the radio-lobes require (this 
estimate, among other things, depends on the proton energy den- 
sity, still a very poorly known quantity). Another inference comes 
from blazars emitting X-rays at large (~ 10-100 kpc) distances as 
observed by Chandra. For them the leading emission model (e.g. 
Tavecchio et al. 2004, 2007; but see. e.g., Kataoka et al. 2008) re- 
quires that the jet is still relativistic at those scales (with F-factors 
similar to the ones derived in the inner regions) and this in turn 
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Figure 4. From top to bottom: distributions of the jet power in the form of 
bulk kinetic power of cold protons (Pp, assuming one proton per emitting 
electron); emitting electrons (Pc, including their average random energy 
(7}moC'^); Poynting flux (Pb), radiation (Pr). The bottom panel shows the 
distribution of the luminosities of the accretion disk. The shaded area in 
all panels corresponds to sources with only upper limits on L^j . 



suggests that the jet has not lost much of its power in producing 
radiation. 

Consider now BL Lacs: we still have that Pr ~ <; Pb, but 
now also Pp is of the same order. This means that we are using 
virtually all the available jet power to produce the radiation we see. 
This is also illustrated in Fig.[5] Thus for BL Lacs we either assume 
that not all electrons are accelerated, allowing for an extra reservoir 
of power in bulk motion of the protons, or, more intriguingly, we 
conclude that the jet noticeably decelerates. The latter option is in 
agreement with recent findings on the absence of fast superluminal 
motion in TeV BL Lacs (Finer & Edward 2004; Finer, Fant & Ed- 
wards 2008), with the absence of strong extended radio structures, 
and with the result of the Chandra observations of extended X-ray 
jets at large scales, showing sub-luminal speed at large scales (e.g. 
Worrall et al. 2001; Kataoka & Stawarz 2005). Moreover, this is- 
sue of jet deceleration of BL Lac jets has been debated recently 
on the theoretical point of view (Georganopoulos & Kazanas 2003; 
Ghisellini, Tavecchio & Chiaberge 2005). 

We conclude that the jet of FSRQs are powerful, matter dom- 
inated and transforming a few per cent of their kinetic power into 
radiation. The jet of BL Lacs are less powerful, with the different 
forms of power in rough equipartition, transforming a larger frac- 
tion of their kinetic power into radiation, and probably decelerating. 
Despite these different characteristics, there is no discontinuity be- 
tween FSRQs and BL Lacs. All the different properties can be ex- 
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tons (fc, mid panel) and in magnetic fields (eg, bottom panel) as functions 
of Pjot = Pp + Pc + Pb - Different symbols refer to different 7-ray lu- 
minosity bins, as in Fig|6] For comparison, we show (little grey circles) the 
blazars studied in Celotti & Ghisellini (2008). 



plained with the difference in jet power accompanied by a different 
environment, in turn caused by a different regime of accretion. This 
important point is discussed below. 

7.1 Jet power vs accretion luminosity 

The availability of the SwiftfU\(JT data for many of our blazars 
made possible to estimate the accretion disk luminosity for sev- 
eral of them. We can then discuss one of the crucial problem in jet 
physics: the disk/jet connection. 

Fig.|6]shows what we think is the main result of our work (see 
also Maraschi & Tavecchio 2003, Sambruna et al. 2006 for earlier 
results). The top panel shows Pr as a function of the accretion disk 
luminosity Ld, while the bottom panel shows Pjct vs Ld. The dif- 
ferent symbols correspond to different bins of the observed 7-ray 
luminosity L-,, as labelled. 

Consider first the top panel: 

• Blazars with different form a sequence in the Pr-L^ plane. 
That L-y correlates with Pr is not a surprise, since we already knew 
(from EGRET) that the 7-ray luminosity is dominating the bolo- 
metric output. What is interesting is that the most luminous 7-ray 
blazars have also a more powerful accretion disk. 

• For our analysis we have considered the average value of the 
7-ray luminosity during the 3 months survey. Then it should not 
be extreme, given the large amplitude and rapid variability shown 
by blazars, especially at high energies. In other words, the shown 
Pr is more indicative of an "average" state, not of an extremely 
high state, even if, in a flux limited sample of variable sources, like 
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the Fermi one, sources in high states are always over-represented. 
Variability of Pi is however an issue, and we can consider that the 
single blazar can vary at least by a factor 10-30 around the shown 
Pr. This contributes to the somewhat large scatter around the Pi- 
Ld relation. 

• Considering only FSRQs, we have that a least square fit yields 
log Pr,45 = 0-73 log I/d.45 ^ 0.36, with a probability for the cor- 
relation to be random of P = 5 x 10~*. The same least square fit 
yields log Ld,45 = 0.65 log Pr,45 + 0.82, indicating that a slope 
around unity is consistent with the data. Since both Pr and Ld de- 
pends on redshift, we have also applied a partial correlation anal- 
ysis, as explained in Padovani (1992). Using Eq. 1 of that paper, 
we have verified that Pr and Ld, once the redshift dependence is 



excluded, still correlate, although the probability to be random in- 
creases to P = 2 X 10"''. 

• Line-less BL Lacs are shown with their corresponding upper 
limits on Ld. These are nevertheless important, showing that they 
must deviate from the general trend defined by FSRQs (see also 
Maraschi & Tavecchio 2003). 

Consider now the bottom panel, showing Pjot vs Ld. As men- 
tioned, the jet power of FSRQs is dominated by the bulk motion of 
cold protons, while in BL Lacs it is more equally distributed among 
electrons, protons and magnetic field. Also in this plane the more 
7-ray luminous blazars have the most powerful jet. 

Considering only FSRQs, we have that a least square fit yields 
log Pjct,45 = 0.62 log Ld,45 + 1.07, with a probability for the cor- 
relation to be random of P = 6 x 10^^. The same least square fit 
yields log Ld,45 = 0.56 log Pjct,45 + 0.14, indicating that a slope 
around unity is consistent with the data. Excluding the dependence 
of redshift by applying a partial correlation analysis, the probability 
to be random increases to P = 3.4 x lO"'^. 

Also in this plane the line-less BL Lacs (with upper limits 
for Ld) deviate from the trend defined by FSRQs, implying that 
their jet is much more powerful than the luminosity emitted by their 
accretion disks. 

We now assume, as an ansatz, that the jet power is always of 
the order of ALnC^, for FSRQ as well as for BL Lacs (see Ghisellini 
& Tavecchio 2008 for a discussion). This allows to estimate the 
accretion rate for BL Lacs independently from their (invisible) disk 
luminosity. Furthermore we can form the ratio Min/MBdd given 
by 



Min 

MEdd 



1.3 X 1038 (M/Mq) 



(5) 



For FSRQs we have used Min and M derived from our modelling 
(Afin is given by ALn = Ld/(r}(?), with r) = 0.08 for all sources). 
For BL Lacs we simply set Mh-i ~ P)ot/(?- The resulting distri- 
butions are shown in Fig. |7] Since for BL Lacs the mass we have 
used is very uncertain, although in agreement with other indepen- 
dent estimates, we show, beside the values of ALn/A/Edd obtained 
using the masses listed in Tab.|4](thick solid line), also the distribu- 
tion obtained by adopting the same mass of 10* (shaded cyan) and 
IQ^Mq (shaded gray) for all BL Lacs. 

Fig.|7]shows that there is a "divide" between BL Lacs and FS- 
RQs occurring at Afin/ALEdd ~ 0.1, equivalent to Ld/LEdd ~ 
0.01, in striking agreement with the value proposed by Ghisellini, 
Maraschi & Tavecchio (2009) and very similar to the value pro- 
posed by Ghisellini & Celotti (2001) for the division between FR 1 
and FR 2 radio-galaxies, based on completely different arguments 
(see also Xu, Cao & Wu 2009). This division can be readily inter- 
preted as the change in the accretion regime of the disk, becoming 
radiatively inefficient when ALn is less than ~ 10% of the Ed- 
dington value (and Ld is less than ~ 1% of LEdd)- Notably, there 
are hints that similar results holds for radio-quiet sources (e.g. Ho 
2009; see also Ho 2008 for review). 



8 SUMMARY OF IfESULTS 

We have studied the entire sample of blazars detected during the 
first 3-months survey of Fermi and have a known redshift and a 
reasonable data coverage of their SED. By studying the resulting 
85 objects we have found the following main results: 

• The simultaneous or quasi-simultaneous Swift observations 
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tom). The assumed black hole mass is listed in Tab. |4] For BL Lacs, given 
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histograms, respectively) for all the sources. 



for a large fraction of our sources allowed to have an unprece- 
dented view on the optical to 7-ray SED of blazars. In addition, the 
optical-UV data were very important to separate the thermal emis- 
sion produced by the accretion disk from the beamed non-thermal 
continuum. In this way, for FSRQs, we could estimate the black 
hole mass and the accretion rate. This in turn allowed to study the 
connection between the power of the jet and the luminosity emitted 
by the accretion disk. We found that they correlate. 

• The estimated black hole masses are in the range between 10* 
and several times 10^ solar masses for FSRQs. For BL Lacs the 
poorly constrained masses are in the range 10* — 10^A/q. These 
values are consistent with those found in the literature for the same 
objects, but existing estimates vary. 

• The luminosity emitted by the accretion disks of FSRQs is 
between I and 60 per cent of the Eddington one. Upper limits to 
the disk emission of BL Lacs indicate Ld/iEdd < 10"^. 

• The "divide" between FSRQs and BL Lacs, in terms of the ac- 
creting mass rate (Ghisellini, Maraschi & Tavecchio 2009), is fully 
confirmed. It occurs when the accretion mass rate becomes smaller 
than 10 per cent of the Eddington one, or, equivalently, when the 
disk luminosity becomes smaller than 1 per cent Eddington. 

• The 7-ray luminosity is a good tracer both of the accretion 
disk luminosity (for FSRQs) and of the jet power (for all blazars). 

• As for the jet emission processes, the EC component almost 
always dominates the emission beyond the X-ray band in FSRQs, 
with the SSC contributing to soft and mid-energy X-rays in some 
cases. In BL Lacs, most of the sources can be fitted by a pure SSC 
model, but some of them require an extra component when the sep- 
aration, in energy, of the synchrotron and Compton peaks is too 
large. This can be provided by a spine/layer structure of the jet, that 
avoids the need of extremely large F-factors. 

• The seed photons for the EC mechanism can be provided by 
a fairly standard BLR, as assumed here and in the "canonical" sce- 



nario for powerful blazars. The majority of FSRQs dissipate within 
the BLR, while 4 of them are better fitted assuming a dissipation re- 
gion between the BLR and an IR emitting torus, at distances greater 
than the BLR. 

• The jet dissipation region is located between a few hundred 
and a thousand Schwarzschild radii for all sources. 

• Bulk Lorentz factors are in the range 10-15. 

• The magnetic field in the emitting region of FSRQs is between 
1 and 10 G, and 10 times less for BL Lacs, on average. 

• Jets in FSRQs must be matter dominated, while in BL Lacs 
there can be equipartition between the power in bulk motion of the 
emitting electrons, cold protons and magnetic field. 

• In FSRQs it is likely that the electron-positron pair component 
is negligible. If so, the jet power in these sources is dominated (by 
a factor 10-30) by the cold proton component, and it is a factor 
^3-5 larger than the luminosity emitted by the accretion disk. The 
outflowing mass rate is around a few per cent of the accreting mass 
rate. 

• We confirm that the SEDs of blazar form a sequence, ex- 
plained in terms of different radiative cooling suffered by the elec- 
trons, with higher energy electrons present in jets of lower power. 



9 DISCUSSION 

The results listed in the previous section confirm earlier findings, 
largely based on blazars with an EGRET detection, and/or detec- 
tions in the TeV band. Because of the factor ~ 20 better sensitivity 
of FennifLKT with respect to EGRET, we are now starting to ex- 
plore sources that are not in "extraordinary" bright states in 7-ray 
band. Sources in our sample should be closer to the average state 
of blazars, even if, given the still limited sensitivity and the large 
amplitude variability (even by factor 30-100) our blazars are most 
likely emitting above their average. 

These Fermi blazars confirm that the jet of blazars form a se- 
quence whose main parameter is their emitted luminosity. This may 
seem strange, given the strong dependence of the observed lumi- 
nosity to the Doppler beaming, then on the viewing angle. On the 
other hand, the results of our model fitting show that our blazars 
are all viewed at small angles, with no misaligned jet entering the 
sample. Misaligned sources, therefore, are fainter than the current 
Fermi blazars and should appear in deeper catalogues. 



9.1 Black hole mass and the blazar sequence 

There is, in our opinion (Ghisellini & Tavecchio 2008), another im- 
portant parameter, besides the jet power, controlling the look of the 
emitted SED and its bolometric luminosity: the mass of the black 
hole. Blazars with small black hole masses and accreting at a rate 
greater than a critical value should be "red" (i.e. they should have 
relatively small peak frequencies and large Compton dominance) 
even if their observed bolometric luminosity is relatively small, 
contrary to what the simplest version of the blazar sequence would 
predict. Fig. [8] shows the observed bolometric luminosity (as de- 
rived by the model) as a function of the black hole mass estimated 
in this paper. Empty circles are FSRQs with estimated black hole 
masses and accretion luminosities, filled circles are BL Lacs with 
only an upper limit on their disk luminosities, and whose black hole 
mass is uncertain. We also show the range of black hole masses ex- 
isting in the literature and reported in Tab. |6] Within FSRQs there 
is indeed a tendency for larger luminosities to correspond to larger 
black hole masses. Viceversa, below Lobs = 10*^ erg s^^ all black 
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Figure 8. The observed bolometric luminosity produced by the jet (calcu- 
lated from the fitting model) as a function of the black hole mass. Empty 
circles are FSRQs with estimated black hole masses and accretion lumi- 
nosities, filled circles are BL Lacs with only an upper limit on their disk 
luminosities, and whose black hole mass is uncertain. We also show the 
range of the estimates of the black hole mass found in the Uterature. 



hole masses are smaller than lO^M©. Fig.[8]shows that there can be 
"red" blazars with Lobs similar to (bluer) BL Lacs, but this happens 
when their black hole mass is relatively small. 

Having the distribution of all relevant physical parameters, we 
can do the exercise to construct the "average" SED of FSRQs and 
BL Lacs, respectively, of our sample. This is illustrated by Fig.|9] 
for which we have used for the average FSRQs and BL Lacs the 
parameters listed at the end of Tab.|4] Note that in our sample there 
are no "extreme" TeV BL Lacs, since, as discussed in T09, these 
sources have so large Compton peak frequencies to make difficult 
a detection by Fermi. So Fig.|9]corresponds to the average BL Lac 
detected by Fermi, and not to the average BL Lac in general. 



9.2 Jet power and accretion 

The jet powers derived here are large, reaching, at the high end 
of their distributions, values greater than 10 times the disk lumi- 
nosity (see the upper panel of Fig. llOb . For these extreme objects 
Pjet>MinC'^. This is admittedly a model-dependent statement. It 
assumes a one-zone leptonic model and that there is one proton 
per emitting electron. A more robust statement is that the amount 
of power the jet spends to produce and carry the non-thermal ra- 
diation is also very large, being in some case equal to the disk lu- 
minosity produced by accretion (see the bottom panel of Fig. 1 10b. 
and more often a factor ~3-10 smaller. Consider that the jet power 
cannot be simply be represented by Pri if this were the case, the 
entire jet power is used to produce the radiation we see, so that the 
jet would decelerate significantly. Instead, the jet must continue to 
be relativistic up to large distances, as required by the existence 
of strong radio lobes and the X-ray radiation seen by the Chan- 
dra satellite at distances (from the black hole) of hundreds of kpc. 
Therefore a reasonable lower limit on Pjct should be a factor 3-10 
greater than P^. 
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Figure 9. The average SED of FSRQs and of BL Lacs without any sign of 
disc emission. These SED have been constructed taking the (logarithmic) 
average of the parameters of the sources belonging to the two subclasses 
(see Table|4). The top and the bottom panels show the fluxes and luminosi- 
ties, respectively. The shown frequencies are calculated in the rest frame of 
the source for the luminosity plot, and are the observed ones for the flux 
plot. The assumed redshift are 2 = 1 for (FSRQ) and z = 0.1 for (BL 
Lac). 
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Figure 10. Top panel: the distribution of the ratio between the total jet 
power Pjct and the accretion disk luminosity L^. Bottom panel: The dis- 
tribution of the ratio between Pr (the power carried by the jet in the form 
of produced radiation) and the accretion disk luminosity . Both distribu- 
tions are for FSRQs only. 
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What is then the source of the power of the jet? Is it only the 
gravitational energy of the accreting matter or do we necessarily 
need also the rotational energy of a spinning black hole? We here 
discuss two possible alternatives, that can both explain our results, 
but are drastically different for the ultimate energy source for the 
jet. 

Jets powered by accretion only — Jolley et al. (2009), building up 
on previous ideas put forward in Jolley & Kuncic (2008), propose 
that, in jetted sources, a sizeable fraction of the accretion power 
goes to power the jet. As a result, the remaining power for the disk 
luminosity is less than usually estimated by setting Ld ~ jyMinC^, 
with rj ~ 0.08-0.1. This implies that the mass accretion rate 
needed to sustain a given Ld is larger than what we have estimated. 
Also the total accretion power is larger, and it is sufficient to explain 
the large jet power we have derived. Assume that the total power 
extractable from the accretion process is r;tot A/inC^, and that a frac- 
tion r]d (Vj) of MinC^ is used to produce the disk luminosity (the jet 
power). We have: 



velocity, we have 



VtotMinC 

rjtot 



rjdMinC 



(6) 



Our results imply i) ryj > rjd and ii) that rjj/rid ~ const in different 
blazars to account for the observed Pj-Ld correlation. 

Jets powered by the black hole spin — The rotational energy of a 
maximally spinning black hole is 29% of the hole rest mass energy 
(i.e. up to 5 X lO^^A/g erg), amply sufficient to power a strong 
jet for its entire lifetime. In principle, in this case one can have 
Pj > -Ld, given a sufficiently efficient way to extract the energy of 
the spinning black hole. In this case we can "decouple" Pj and Ld, 
since they have a different energy source. 

On the other hand, for FSRQs, we do see a relation between 
the Pj and Ld, and at first glance this seems to suggest that it is the 
accretion, not the spin, to power the jet. We can envisage a possible 
solution to this apparently contradictory issue, by linking the ex- 
traction of the hole rotational energy to the accretion process. The 
main idea is the following: the energy density pov^ of the accret- 
ing material close to the black hole horizon can sustain a maximum 
magnetic energy density Bq / (Svr) of the same order (see also Ghis- 
ellini Sl Celotti 2002). Here v^, is the circular (keplerian) velocity 
of the matter. The magnetic field sustained by the accreting matter 
can then tap the rotational energy of the hole. The mechanism able 
to do this task is the Blandford-Znajek (1977) process (hereafter 
BZ), whose efficiency has been debated in recent years (e.g. Mod- 
erski & Sikora 1996; Ghosh & Abramowicz 1997; Livio, Ogilvie 
& Pringle 1999; McKinney 2005; Garofalo 2009; Krolik & Haw- 
ley 2002). Here, for simplicity, without entering in the technical 
discussion on the efficiency of this mechanism, we assume that the 
jet power is of the form: 



; 2 n2 Bo 

Ka nHo c 

Stt 



(7) 



where a ^ 1 is the dimension-less angular momentum of the hole, 
Ro is some fiducial distance of the order of the black hole horizon, 
and Bo is the magnetic field at that radius. Eq. [7] is nothing else 
than a Poynting flux. The factor k includes our uncertainties about 
the efficiency of the BZ process. 

Assume that at a distance Rq from the black hole, the disk has 
a height Ho above the equatorial plane. If Vi is the radial in-fall 



Min = ^TvRoHopovn 



(8) 



AuRqHovb^ 

If the energy density of the magnetic field is a fraction ea of the 
kinetic energy density of the matter orbiting around the black hole, 
we have 



Bo 1 2 

8^ = 2^^^°"^ 



(9) 



If a magnetic field with the same magnitude is threading the spin- 
ning hole, then we have 



8(Ho/Ro) I3r 



8v{Ho/Ro) ^ 



(10) 



Close to the gravitational or Schwarzschild radius, we may set 
Ho/Ro<l. The ratio P^/Pb. can be slightly larger than unity, de- 
pending on viscosity. For ee ~ 1, the jet power is maximum, and 
for k not much less than unity can be of the same order of the disk 
luminosity. The (rather strong) requirement that our data are pos- 
ing on the BZ mechanism is therefore on its efficiency, that must be 
large. 
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Log u [Hz] 

Figure 11. SED of PKS 0166-219, 0133+47 (DA 55) and PKS 0142-278, 
together with the fitting models, with parameters listed in Tab.|4] Fermi and 
Swift data are indicated by dark grey symbols (red in the electronic version), 
while archival data (from NED) are in light grey. The short-dashed line is 
the emission from the IR torus, the accretion disk and its X-ray corona; the 
long-dashed line is the SSC contribution and the dot-dashed line is the EC 
emission. The solid light grey line (green in the electronic version) is the 
non thermal flux produced by the jet, the solid dark grey line (blue in the 
electronic version) is the sum of the non-thermal and thermal components. 
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Figure 12. SED of PKS 0208-512, B2 0218+35, PKS 0332^03 and PKS 
0537^41. Symbols and lines as in Fig. 1111 
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Figure 13. SED of B3 0650+453, TXS 0716+332 and 0814+425 (OJ 425). 
For the latter source we show a pure SSC model and an EC one, as indicated. 
Symbols and lines as in Fig. 1111 
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Figure 14. SED of 0954+556 (4C 55.17), S4 1030+61. PKS 1055+018 and j-jg^^e 15. SED of PKS 1127-145, PKS 1144-379. 1156+295 (4C 29.45) 

PKS 1057-79. Symbols and lines as in ¥ig.\T^ and 1226+023 (3C 273). Symbols and lines as in Fig.O 
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Figure 16. SED of PKS 1244-255, 1253-055 (3C 279), B2 1308+32 and 
PKS 1352-104. Symbols and lines as in Fig.fTT] 
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Figure 17. SED of PKS 1508-055, PKS 1510-089, PKS 1622-253 and S5 
1803+784. Symbols and lines as in Fig.QT] 
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Figure 18. SED of TXS 1846+322, S4 1849+67, PKS 1908-201 and TXS 
1920-21 1. Symbols and lines as in Fig. II II 



Figure 19. SED of 2141+175 (OX 169), PKS 2144+092, B2 2155+31 and 
PKS 2201+171. Symbols and lines as in Fig.fTT] 
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Figure 20. SED of PKS 2204-54, 2230+114 (CTA102) and PMN 2345- 
1555. Symbols and lines as in Fig. 1111 
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